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ABSTRACT 


Recent evidence on the metal content of the high-redshift Lya forest seen 
in quasar spectra suggests that an early generation of galaxies enriched the 
intergalactic medium (IGM) at z ^ 5. We calculate the number of supernovae 
that need to have taken place to produce the observed metallicity. The 
progenitor stars of the supernovae should have emitted ~ 20 ionizing photons 
for each baryon in the universe, i.e., more than enough to ionize the IGM. We 
calculate that the rate of these supernovae is such that about one of them should 
be observable at any time per square arc minute. Their fluxes are, of course, 
extremely faint: at z = 5, the peak magnitude should be A = 27 with a duration 
of ~ 1 year. However, these supernovae should still be the brightest objects in 
the universe beyond some redshift, because the earliest galaxies should form 
before quasars and they should have very low mass, so their luminosities should 
be much lower than that of a supernova. 

We also show that, under the assumption of a standard initial mass function, 
a signihcant fraction of the stars in the Galactic halo should have formed in 
the early galaxies that reionized and enriched the IGM, and which later must 
have merged with our Galaxy. These stars should have a more extended radial 
distribution than the observed halo stars. 

Subject headings: Galaxy: halo - galaxies: formation - large-scale structure of 
universe - quasars: absorption lines - supernovae: general 
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1. Introduction 

Ever since the discovery of the first high-redshift quasar (Schmidt 1965), quasars have 
maintained their title as the objects with the highest known redshift; the present record 
holder is a quasar at z = 4.89 (Schmidt, Schneider, & Gunn 1991). Nevertheless, the 
highest known redshifts of galaxies have followed closely behind, with bright radio galaxies 
having been found up to 2 ; = 4.45 (Lacy et ah 1994; Rawlings et al. 1996); more recently, 
galaxies with high star formation rates have been identihed from interstellar absorption 
lines at 2 : = 2 to 3 (Steidel et al. 1996) and from the Lyo; emission line at z = 4.55 (Hu & 
McMahon 1996). 

In fact, if quasars are related to supermassive black holes that formed in the centers of 
high-redshift galaxies, we should expect that many galaxies already existed before the hrst 
quasars appeared. In any ‘bottom-up’ theory where the observed structure in the universe 
forms by hierarchical gravitational collapse, and the primordial density fluctuations extend 
to sufficiently small scales, the first galaxies to form must have had much smaller masses 
than the present galaxies. The first stars should have formed in systems with velocity 
dispersions of ~ lOKms”^ or lower, corresponding to the lowest temperatures (T ~ 10^ 

K) that allow cooling and dissipation of the gas by atomic processes (systems with even 
lower virial temperatures can cool and dissipate through molecular hydrogen, but this 
cooling process should be suppressed by photodissociation of the molecules after emission 
of a number of UV photons that is much smaller than that needed to reionize the universe; 
see Haiman, Rees, & Loeb 1996). These systems would be very unlikely to form quasars, 
because even a small fraction of their baryons turning into stars should provide sufficient 
energy (via ionization, stellar winds or supernovae) to expel the remaining gas from the 
shallow potential well (e.g., Couchman & Rees 1986; Dekel & Silk 1986). Deeper potential 
wells, forming at later epochs, are probably needed to form supermassive black holes in 
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galactic centers. 

Even if all the baryons were converted into stars very efficiently in these early dwarf 
galaxies, with a baryonic mass Mb ^ IO^Mq, their total stellar luminosity would be much 
smaller than in L* galaxies at present, simply due to their small mass. Since only a small 
fraction of the baryons in these systems is likely to turn into stars before the gas is ejected, 
the total stellar mass in the first galaxies to form in the universe should be much smaller 
than IO^Mq. This implies that a supernova in one of these first galaxies to form will be far 
brighter than the galaxy itself. Thus, the brightest probes of the era when the reionization 
of the intergalactic medium (IGM) started should be supernovae in very small galaxies, 
caused by the death of the same stars responsible for the hrst ionizing photons. In this 
paper, we shall estimate the number of supernovae that should have taken place in these 
galaxies and should be observable at very high redshift, and their apparent magnitudes. 


2. Rate of High Redshift Supernovae 

The scenario where these first small galaxies caused the reionization of the universe 
is strongly supported by recent evidence that the metal abundance in the Lya forest 
absorption lines with N^i ^ 10^"^ cm“^ is typically Z r'U 10-2^0, from observations of 
CIV lines (see Tytler et ah 1995, Songaila & Cowie 1996 and references therein). There 
is relatively little uncertainty in the number of UV photons that were emitted by the 
stars that produced a given mass of heavy elements, because the heavy elements originate 
from the supernovae resulting from the same stars that emit most of the ionizing photons 
(although the C/0 ratio is more uncertain because carbon is more abundantly produced 
in lower mass stars). According to the most recent calculations (Madau & Shull 1996 
and references therein), the energy of Lyman continuum photons emitted is 0.2% of the 
rest-mass energy of the heavy elements produced. Thus, the energy emitted in ionizing 
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photons per baryon is 0.002 = 2Z MeV, where Z is the average metallicity of all 

baryons in the nniverse, so only Z = 10“^ is needed to have emitted one ionizing photon for 
each baryon. \i Z = IO^^Zq = 2 x 10“*^, then 20 ionizing photons must have been emitted 
per baryon when the heavy elements were made. Furthermore, if these were the photons 
responsible for reionizing the universe, then each baryon must have recombined 20 times on 
average during the reionization epoch. This is a reasonable number, because a fraction of 
these photons were probably absorbed in the systems where the stars were born before the 
gas was expelled, and those that escaped could also have been absorbed in Lyman limit 
systems before the universe became transparent. Thus, there is no need to invoke ejection 
of gas by more massive galaxies that can accrete the ionized IGM to explain a metallicity 
Z = 10-2^0. 

We can now calculate the number of supernovae that were required for enriching the 
gas in the IGM to the average metallicity of Z = 2 x 10“"^ that is observed in the Lya 
forest at X ~ 3. This number should depend only on the IGM density and the supernova 
yields, and should be independent of any other details related to the theory for galaxy 
formation and the type of galaxies that ejected the enriched gas. Recent simulations of 
cold dark matter models show that the absorption lines in the Lya forest can be identihed 
with the IGM, with density fluctuations caused by gravitational collapse, and that most of 
the baryons should be in the IGM in these models (Gen et al. 1994, Hernquist et al. 1996, 
Miralda-Escude et al. 1996). Thus, it is reasonable to assume that the high-z supernovae 
enriched all the baryons in the universe to a mean at metallicity at least as high as that of 
the Lya forest. 

Since each supernova produces an average of ~ IMq of heavy elements (with 
uncertainties depending on the assumed initial mass function and supernova models; 
see Woosley & Weaver 1995), this implies that a supernova took place at high redshift 
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for each SOOOMq of baryons in the universe. We shall assume a high baryon density 
Qb = 0.025h~^, in agreement with the primordial deuterium abundance measured by 
Buries & Tytler (1996) and by Tytler, Fan, & Buries (1996). Notice that this implies 
that most of the baryons at the present time are dark, so many more baryons than 
those we observe in galaxies had to be enriched at high redshift. Assuming the hi = 1 
cosmological model, the total mass of baryons in a redshift shell of width around us 
is Mb = {6c^ilb)/{GH q) [1 — (1 + /{\ + where hfo is the present Hubble 

constant. With the above rate of supernovae per baryon mass (assumed to take place 
within the epoch corresponding to the redshift shell Az), and taking into account that the 
supernovae within the shell would be seen by us over a time interval lS.zj{\ + z'f'l‘^ ^ the 
total supernova rate observed over all the sky is 

Rsuv = 1-8 X 10®h"^ [1 - (1 + z)~^^‘^f yr"^ . (1) 

or, for 2 : ~ 5, about 1 supernova per square arc minute per year. 


3. Apparent Magnitudes 

Most of these supernovae should be Type II which, if the progenitor is a red supergiant, 
have a plateau of the luminosity in their lightucurves from 1 to 80 days after the explosion, 
with L ~ 3 X 10^^ erg/s (Woosley & Weaver 1986). A note of caution should be made here, 
in that the low-metallicity progenitors of these early supernovae could be very different from 
the high metallicity counterparts. As illustrated by the case of SN1987A, it is probably not 
possible at this stage to predict the type of supernovae we should expect from this first 
generation of stars; in particular, the possibility that some supernovae might reach higher 
intrinsic luminosities than regular Type IIs should be kept in mind. A duration of 80 days 
would be redshifted to more than a year at redshifts z > A, where the supernovae from 
stars responsible for the reionization of the universe should occur. At any random time 
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there should therefore be one or more supernovae per square arc minute visible in the sky. 
Of course, many more supernovae should have occurred in more massive galaxies at later 
epochs (producing the metals in stars and interstellar gas at present), and possibly in small 
systems that ejected their gas and continued to enrich the IGM. 

The main difficulty in detecting these supernovae is obviously the extremely faint 
flux expected. Supernovae should be the brightest objects in the universe at very high 
redshift, but they are of course much fainter than quasars and, as the redshift increases, the 
bolometric fluxes decrease at least as rapidly as {1 + z)"^. One should point out, however, 
that when observing at a hxed, long wavelenth such that the supernovae are observed on 
the Rayleigh-Jeans part of the spectrum, the flux actually becomes brighter as 1 + ; 2 , in the 
limit of high redshift. 

To estimate the apparent magnitudes of the supernovae, we assume a blackbody 
spectrum, which is a sufficiently close approximation for our purpose. In Figure 1 we 
show the apparent magnitudes in several bands for a luminosity L = 3 x 10^^ erg/s, and 
temperatures T = 25000 K and T = 7000 K. The high temperature is reached rsj 1 day 
after the explosion, when the luminosity drops to the value in the plateau part of the 
lightcurve; during the next few days the temperature cools, reaching a value near 7000 
K after about a week, and then it stays constant until the luminosity starts decreasing. 
Immediately after the explosion, when the shock reaches the surface of the star, the 
luminosity and temperature can be much higher and the apparent flux can be brighter by 
a factor ~ 10 relative to the values in Fig. 1, but this phase only lasts for ~ 1 hour. The 
apparent magnitudes have been calculated for the hi = 1 model with Hq = 70 Km/s/Mpc, 
from the flux at the central wavelengths of the bands, using the central wavelengths and 
zero-magnitude fluxes given in Allen (1973). The high-redshift cutoff in the curves in Figure 
1 indicate the redshift where the supernova light would be absorbed by the Lya forest. 



At high redshift, the supernovae must obviously be searched in the infrared. The 
faintest galaxy surveys from the ground have reached magnitudes A' ~ 25 (Cowie et 
ah 1994). From Fig. 1, type II supernovae would have similar magnitudes at ~ 2.5, 
although searching for variable objects should require more telescope time than simple 
object identihcation. In order to hnd supernovae at redshifts higher than known quasars, 
fainter fluxes by a factor of 5 — 10 need to be detected. This might be achieved with 
implementation of adaptive optics on large ground-based telescopes; in the longer term, the 
New Generation Space Telescope should certainly be capable to observe these supernovae 
(see Mather & Stockman 1996). We emphasize again the large uncertainty in predicting the 
types and absolute magnitude of supernovae from these early generation of stars. 

One possibility to detect these supernovae before more powerful telescopes and cameras 
in the infrared can be built is to use the magnifying power of gravitational lensing in 
galaxy clusters. The deflection angles of the most massive clusters of galaxies are as large 
as 6 ~ 30'', with critical lines of total length of several arc minutes. The cross section for 
magnifying a source by more than a factor A is TTb^jA^, or 0.01 square arc minutes for 
A = 10. Thus, any rich lensing cluster should have a 1% chance of having a high-redshift 
supernova magnihed by a factor larger than 10 at any time. The highly magnihed images 
would always appear in pairs around critical lines, and would be simple to identify only 
from their positions and colors given a lensing model of the cluster (see Miralda-Escude & 
Fort 1993). Later, variability would have to be detected to distinguish the supernova from 
a faint, compact galaxy. 


4. The Present Distribution of Population III Stars 

We have suggested in this paper that, under reasonable assumptions, supernovae should 
be the brightest objects in the universe beyond some redshift, in particular during the 
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early phases of the reionization. The supernovae might therefore be the first observational 
evidence we shall have of this epoch, when the very faint apparent magnitudes expected 
are observable. The other observable signature of this epoch may be the 21cm absorption 
or emission by the neutral intergalactic gas (Scott & Rees 1994; Madau, Meiksin & Rees 
1996). 

The UV and heavy elements abundance inferred from quasar absorption lines allow us, 
as we have seen, to draw quantitative conclusions about the minimum number of high-mass 
stars formed beyond z = 5: to produce a metallicity Z = 2 x 10“"^ requires one supernova for 
each ~ 5 OOOM 0 of baryons. This inference is quite robust, being insensitive to the details of 
structure and galaxy formation at high redshift, which of course depend on the cosmological 
assumptions. However, the total mass of stars formed depends on the IMF, and is therefore 
much more uncertain. For a standard IMF, ~ IOOMq of stars need to be formed to produce 
one supernova, so 2% of the baryons should have been turned into stars by the time the 
IGM reached this level of enrichment. Notice that this is equivalent to 20% of the observed 
baryons in galaxies today, given our adopted value of Vlb which implies that only ~ 10% of 
the baryons are in known stars and gas in galaxies. It is of course quite possible that the 
IMF was different for these early stars, given the different physical environment (a higher 
ambient temperature, absence of heavy elements to act as coolants and provide opacity, 
and no signihcant magnetic fields). Direct clues to the slope of the high-mass IMF may 
come from (for instance) the relative ionization levels of H and He and heavy elements, 
which depend on the background radiation spectrum shortward of the Lyman limit, or from 
relative abundances of heavy elements relative to carbon. Conceivably, all the early stars 
might be of high-mass, so that no coeval low-mass stars survive; at the other extreme, the 
early IMF could have been much steeper than the standard one, in which case there could 
be many pregalactic brown dwarfs. 
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Would any of these “Population III” stars be observable today? Let us consider the 
observational consequences of the simplest assumption: that the early IMF was the same as 
in the solar neighborhood. In that case, most of the present luminosity from the Population 
III would arise from red giants and stars at the tip of the main-sequence, with M ~ O.SMq. 
Where should these stars be today? 

After the hrst galaxies ejected all their gas back to the IGM, the stars that had been 
formed should have remained in orbit near the center of the dark matter halos. The stars 
then behave as collisionless matter as the halos merge with larger objects, until the present 
galaxies are formed. We would therefore expect that these stars would at present be 
distributed approximately like the dark matter in galactic and cluster halos, and in addition 
there should be some surviving galaxies from that epoch which have not merged into much 
larger objects (or have survived in orbit after merging with a large halo, having escaped 
tidal disruption) and still have the Population III stars in their centers. The halos of stars 
formed in this way around galaxies might be somewhat more centrally concentrated than 
the dark matter, if many mergers take place with only a moderate increase of the halo mass 
at each merger (so that dynamical friction is effective after each merger and it can bring 
the stars near the center of the newly formed halo before tidal disruption occurs). In fact, 
particles that start near the centers of halos that merge tend to end up near the center of 
the merger product (e.g., Spergel & Hernquist 1992). 

The known halo stars have a very steep density prohle, p oc and their total mass 
is M ~ IO^Mq (e.g., Morrison 1993). This mass is comparable to the total mass we would 
expect in the halo in the Population III stars, if the total mass of the halo of our Galaxy 
is 5 X IO^^Mq, with a baryon fraction of 10%, and if 2% of the baryons formed Population 
III stars. Therefore, if the stellar mass function in the hrst galaxies was normal, a sizable 
fraction of the halo stars should have originated there (this is not surprising, because it is 
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derived from the assumption that the halo stars created their own metal abundance). It 
seems difficult that the process of dynamical friction alluded to above can result in the steep 
slope of the halo stars. However, the halo density prohle might become shallower at large 
radius (see Hawkins 1983 and Norris & Hawkins 1991 for current observational evidence on 
this possibility), and a second halo population in the outer part of the galaxy (i? ~ 100 
Kpc) might be the remnant of the Population HI. These halo stars could be found in the 
Hubble Deep Field (HDF, Williams et ah 1996). If the stellar mass of this outer halo is 
IO^Mq, there should be ~ 10® stars near the main-sequence turnoff, i.e., we expect a few 
stars in the HDF (with area 4.4 arcmin^); these would have colors / — H ~ 1.5, J ~ 25 at 
distances of 100 Kpc. From Fig. 2 in Flynn, Gould, & Bahcall, we see that there is at least 
one stellar object with these characteristics in the HDF. 

Several other observations may help to test the existence of the Population HI stars. 
An outer stellar halo would also imply a certain number of high-velocity stars near the solar 
neighborhood. A stellar population may be found in the halos of external galaxies, with 
density prohles similar to the dark matter. Sackett et al. (1994) found a luminous halo 
in the galaxy NGC 5907 with M/L = 500 (i.e., about ten times more light than what we 
expect for the Population HI). Planetary nebulae could also be found in nearby halos of 
galaxies or galaxy groups; several of them were reported recently by by Theuns & Warren 
1996) in the Fornax cluster. 

There is also the possibility that the IMF in the early galaxies produced a large number 
of brown dwarfs. In this case, a large fraction of the baryons could have been turned into 
brown dwarfs, and these could be detected in ongoing microlensing experiments towards 
the LMC (see Paczyhski 1996). If the baryon fraction in the universe is 10%, the optical 
depth of these brown dwarfs toward the LMC could be as high as a few times 10“®. 


Finally, we notice that the metallicity distribution of the Population HI stars is difficult 
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to predict. If only a small fraction of the neutral IGM collapsed to galaxies before the 
reionization, then the gas in these galaxies could reached high metallicities and formed stars, 
and the metallicity could have diluted in the IGM when the gas was ejected. At the same 
time, the metal abundance of the IGM after reionization could be highly inhomogeneous, 
so some galaxies formed later could have very low metallicities. Therefore, it is difficult to 
predict even if the average metal abundance of the Population III stars should be higher, 
lower or similar to the more centrally concentrated halo stars, let alone the distribution of 
these metallicities. 


5. Conclusions 

As the observational techniques improve our ability to detect extremely faint sources, 
and higher redshift objects can be searched for to continue unravelling the history of galaxy 
formation, supernovae should become the brightest observable sources. These supernovae 
created the heavy elements that were expelled to the IGM, and their progenitor stars are the 
most likely sources of the photons that reionized the universe. The expected rates of these 
supernovae, calculating under the assumption of a high baryon density {VL^h? = 0.025), 
and an average metal production oi Z = 10 ~^Zq, is as high as 1 supernova per square arc 
minute per year. To detect the supernovae, the flux limits of the faintest sources detectable 
with our telescopes will probably need to be pushed by another ~ 2 magnitudes, although 
the hrst examples might be discovered at brighter fluxes behind clusters of galaxies, using 
the lensing magnihcation. 

Any low-mass stars that were formed in the hrst small galaxies where these supernovae 
took place should be observable today. We have argued that, if the IMF in these galaxies 
was similar to the present one in our galactic disk, the Population III stars are likely to 
account for a large fraction of the stars in our galactic halo, although most of them should 
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be in an as yet undetected outer halo with a shallower density prohle than the known, inner 
stellar halo. 

We thank Len Cowie, Andy Gould and John Norris for stimulating discussions. JM 
acknowledges support by the W. M. Keck Foundation at IAS. 
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Fig. 1.— Apparent magnitude of a supernova as a function of redshift in different bands, as 
labeled in the figure, for two different temperatures. We assume the supernova spectrum is 
a blackbody with the luminosity and temperature indicated. Supernovae Type II generally 
have a plateau of constant luminosity in their lightcurves with the value assumed in this hgure 
lasting for 80 days. A temperature of 25000 K is reached ~ two days after the explosion, 
and a week later the temperature has dropped to near 7000 K, where it stays constant for 
the next two months. 








